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ABSTRACT 

A simple code for use in analyses of gamma radiation effects 
in laminated materials is described in this note. Simple and good 
geometry is assumed so that all multiple collision and scattering 
events are excluded from consideration. Buildup factors may be 
applied to the results of this simplified analysis, when available, 
to arrive at improved estimates of energy deposition and attenuated 
intensity for actual materials. 

The code proper is capable of handling laminates up to six 
layers. However, for laminates of more than six layers, the same 
code may be used to incorporate two additional layers at a time, 
making use of punch- tape outputs from previous computation on 
all preceding layers. 

Spectrum of attenuated radiation may be obtained as both 
printed output and punch tape output as desired. This code is 
written in "Extended Basic" language, compatible with the GE 26 5 
Time Sharing system. 


I, Theoretical Fundamentals 

Energy deposition resulting from inter- 
action between photons (X-rays, gamma rays) 
and matter has been discussed in some detail 
in Reference (1). The basic results, obtained 
for simple and good geometry in a one-dimen- 
sional representation, are 


Attenuated spectrum at depth x: 
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Energy deposited at depth x: 
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where x 
dm 


= depth into the material, in cm. 
= p dx = differential mass depth 
in gm/ sq cm. 


P 

E 

I 

o 


1(E) 


cr(E), a (E) 
a 


density of material, in gm/cc. 
energy level, in KeV 
fluence, in cal/sq cm, 
incident on the surface of the 
first layer 

spectrum of the incident 
radiation, function of E, 
normalized to unity for E 
ranging from 0 to infinity 
mass attenuation and mass 
absorption coefficients 


To circumvent the need for a detailed defin- 
ition of incident energy spectrum, Planckians 
expressed in terms of equivalent blackbody 
temperatures, in KeV units, are used as a 
convenient means to segregate engineering 
estimates of shielding and thermomechanical 
effects completely from specific identification 
of spectral contents of the primary source, 
i.e. , bomb types, etc. Therefore, the 
quantity 1(E) may be thought of as 

1(E) = I. . P(E/E.) 
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where I. = fluence weighting factor for 

the i-th equivalent black body, 
with Si I. = 1, and 

P($) - normalized Planckian, as 
shown in Figure 1. 

Since I Q appears in both expressions displayed 
above, no loss of generality will result if I Q 
is set equal to unity, corresponding to a 
fluence of 1 calorie per square centimeter. 
Thus a normalized dose Q(x) may be defined 
as 

Q(x) = q(x)/I 

o 

which is a unique property of the material in 
question when subjected to an incident spec- 
trum associated with a specific equivalent 
blackbody temperature, or with a particular 
combination of spectra associated with a 
number of black bodies of different equivalent 
temperatures. The attenuated spectrum is of 
course given by the ratio ^(EJ/E, or simply 
I X (E), since I 0 = 1. 

When materials of various thicknesses 
are placed one after another, either as close- 
ly spaced laminates or as a series of material 
layers with air spaces in between, the above 
computation process can be applied to each 
layer in succession. The attenuated spectrum 
from the first layer becomes the incident 
spectrum for the next layer, and so on. Thus, 
the entire series of computation may be made 
by means of a simple loop; the number of 
layers that may be accommodated becomes 
strictly a question of the size of machine 
storage available. 

The mass absorption coefficients and the 
mass attenuation coefficients are both com- 
plicated functions of energy level, which have 
been tabulated by several agencies. For 
definiteness, the following source data are 
used -- these data are not always consistent 
with each other, but they represent the best 
information available at this time. Logarith- 
mic interpolation is used throughout for 
establishing pertinent values associated with 
each energy level of concern. 

(1) Photon Cross Sections from 0. 001 to 
100 MeV for Elements 1 through 100, 
Los Alamos Scientific Laboratory 
Report LA-3753, November 15, 1967. 

(2) Analytical Approximations for X-ray 
Cros s Sections, Sandia Laboratory, 
Report SC -RR- 66-452, February 
1967. 

(3) X-ray Attenuation and Absorption 
Coefficients, The Boeing Company, 
Repo rt D2 - 1 25 065 - 1 , September 
1966 . 


A new edition of item (1) is in press, and 
will be scanned for necessary updating when 
published. Values to be used can be com- 
posites based on all three of the above ref- 
erences, a sample of which is shown on the 
next page, for the element copper. 

For compounds, both coefficients may 
be approximated by proportional weighting 
in accordance with their atomic abundancy. 
Since the mass coefficients have already been 
referred to unit mass densities, this implies 
that 

T = Z* A x. N i ^X./Z* A x. N i 
1 l / 1 

where Axj is the atomic weight of the ele- 
ment Xj. The molecular formula of the 
compound in question is of course 

TT (x i'N. 

i 

For heterogeneous mixtures of several com- 
pounds, the attenuation problem involved is 
far more complex than can be correctly 
handled by the present simplified approach. 
However, for an engineering estimate, it 
will be assumed that they behave like com- 
pounds with the same proportion of various 
elements by weight. 

II. Description of the Code ri ENEDEP" 

The code written in the Extended Basic 
language for use on the GE-265 Time Shar- 
ing system is named ENEDEP (ENErgy 
DE Position). It is capable of handling lami- 
nations up to six layers of different materials. 
Various combinations of source temperature 
(equivalent blackbody temperature) and inci- 
dent fluence can be processed in series for 
the same material configuration. The output 
of this code provides: 

(1) Identification of layer material by 
name, its density and thickness. 

(2) Depth, mass depth, energy deposi- 
tion per unit mass and local radia- 
tion intensity in tabular form. 

(3) Attenuated spectrum at the back 
side of the last material layer, in 
tabular form and/or in punch tape 
form as desired to be selected by 
appropriate control indices. 


Input data requirement, 
possible modification to accomodate mater- 
ials with a large number of edges in their 
mass absorption coefficients, and a sample 
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problem are discussed in the following 
sections. 

III. Input Data 

Input data for executing the ENEDEP 
code must be presented in the following 
format: 

1200 DATA N, M 

N index for printout and punch tape 
input selection: 

0 no attenuated spectrum informa- 
tion at all 

1 tabulated attenuated spectrum, 
punched tape output optional 

2 punched tape input and tabulated 
attenuated spectrum with optional 
punched tape output 

M integer, denoting the number of 
layers, not greater than 6. 

2000* DATA X(I) N(I), : 

X(I) Thickness of the I-th layer, in 
mils 

N(I) Number of printout desired in the 
I-th layer, integer 

*2000 is a typical line number, which 
may be any number greater than 
1200, followed by tape input 
containing 

2010* "AAAAAAAAAA", W(I), J(I), : 

n AAA AAAAAAA ' 1 alpha-numeric 
string of ten characters, being 
the name of the material of the 
I-th layer 

W(I) Density of the material in the 
I-th layer, in gm/cc 

J(I) Number of sets of mass atten- 
uation and mass absorption co- 
efficients, integer, not greater 
than 35 

*2010 is a typical line number; for 

materials with no edges between 
0. 5 KeV and 200 KeV, there will 
be 22 sets of coefficients on this 
prepared tape; for materials with 
edges, add two sets for each edge 

9900 DATA Sy P y S 2 , P 2> 

S. Incident flu ences, in cal/sq cm 
1 (equal to 1 for most cases) 

P. Equivalent blackbody temperature, 

1 in KeV units 


The usual requirement calls lor five 
sets of S, P combinations, namely, 1,1, 
1,2, 1,4, 1,7, 1,15. However, any 
other combinations pertinent to a par- 
ticular problem can be used. 

On the prepared tape, in addition to what is 
explained above under "201 0", there will be 
other lines of material properties coefficients, 
specified in triplets, giving: 

E(j) Energy coordinate, in KeV 

F(j) Mass attenuation coefficient, a 

G(j ) Mass absorption coefficient, a a 

IV. Materials With More than Six Layers 

When the material under consideration 
has more than six layers, the first six 
layers will be processed as indicated above, 
using an index N = 1 and obtain a punch tape 
output. This tape is then read into a new 
problem, with the following important 
modifications: 

1200 DATA 2, 2 

Only two more layers may be treated, 
due to storage limitation in the GE-265 
Time Sharing system. 

Line 9900 must be modified to read as 
follows: 

9900 DATA S, P 

S, P here must be the same set of 
parameters pertaining to the punch tape 
data being read in as tape input. 

On occasions, more than 35 entries may 
be required to provide an adequate repre- 
sentation of the mass attenuation and mass 
absorption coefficients of one or more 
materials. In such cases, the dimensional 
statement in line 20 must be revised accord- 
ingly, an example of this may be seen in 
section V. Due to limitation in the size of 
data storage capacity, the number of layers 
of materials may also be adjusted. Thus a 
typical modification of line 20 may have the 
following form: 

20 E(4, 42), F(4, 42), G(4, 42) 

V, Sample Problem 
To illustrate the usage of the ENEDEP 


code, a six layer laminate consisting of the 
following materials is analyzed: 

Quartz 

0.013 gm/cm 2 2 mils 

Silicon 

o 

1.07 gm/cm 180 mils 

Quartz 

0.013 gm/cm 2 2 mils 
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Platinum 0. 0022 gm/cm 2 0. 04 mils 

Gold 0. 002 gm/cm 2 0. 04 mils 

Bi 2 Te 3 2,0 gm/cm 2 100 mils 

Due to the high Z -number associated with 
the materials in the last three layers, 42 
entries are required for each of these 
layers. Therefore, the main program is 
modified to handle only four layers, each 
with provision for 42 entries. The main 
program would have to be recycled three 
times to complete the analysis required 
for all six layers. Input data for the first 
four layers used in the first cycle, results 
for the first four layers, input data for the 
fifth layer with tape input of attenuated 
spectrum from the first four layers, out- 
put for the fifty layer, input data for the 
sixth layer with tape input of attenuated 


spectrum from the first five layers, and 
output for the sixth layer are displayed in 
succession at the end of this paper. A 
summary of results for the entire range 
of equivalent blackbody temperature from 
1 to 16 KeV is displayed in Figures 2 and 3. 

VI. Concluding Remarks 

A simple code for determining energy 
deposition dosage and attenuated fluence in 
material laminates induced by gamma ray 
radiation is presented. The average run- 
ning time for one set of fluence and equiva- 
lent blackbody temperature is approximately 
80 seconds (computer time). This code 
therefore provides a very convenient and 
economical tool for parametric studies in 
shielding problems and energy deposition 
analyses. 
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COPPER 


29; A - 63»5h; 


8.9lis 


28 


gi l en in noatin g Point format} powers of 
ten are shown after the letter E. 


Fhton Energy 

(Kev) 


Attenuation Absorption 


1 

2 

3 

1 * 

5 

6 

7 

8 

9 

10 

11 

12 

13 

U* 

15 

16 

17 

18 

19 

20 

21 

22 

23 

21 * 

25 

26 

27 

28 


0.5000 E 00 
0 . 6000 E 00 
0.8000 E 00 
0.9320 E 00 
0.9320 E 00 

0.1000 E 01 
0.1100 E 01 
0.1100 E 01 
0.1500 E 01 
0.2000 E 01 

0.3000 E 01 
0.1*000 E 01 
0.5000 E 01 
0.6000 E 01 
0.8000 E 01 

0.8980 E 01 
0.8980 E 01 
0.1000 E 02 
0.1500 E 02 
0.2000 E 02 

0.3000 E 02 
0.1*000 E 02 
0.5000 E 02 
0.6000 E 02 
0.8000 E 02 

0.1000 E 03 
0.1500 E 03 
0.2000 E 03 


0.1120 E 05 
0.7007 E 01* 
0.3333 E 01* 
O.2237 E 01* 
0.161*6 E 05 

0.1358 E 05 
0.101*6 E 05 
0.1100 E 05 
0.5250 E 01* 
0.21*00 E 01* 

0.7800 E 03 
0.3500 E 03 
0.1900 E 03 
0.1110 E 03 
0.51*00 E 02 

0.1*100 E 02 
0.2960 E 03 
0.221*0 E 03 
0.7600 E 02 
0.3300 E 02 

0.1060 E 02 
0.1*570 E 01 
0.21*10 E 01 
0.11*90 E 01 
0.7080 E 00 

0.1*250 E 00 
0.2080 E 00 
0.11*80 E 00 


0.1120 E 05 
0.7007 E 01* 
O.3333 E 01* 
0.2237 E 01* 
0.161*6 E 05 

0.1358 E 05 
0.101*6 E 05 
0.1100 E 05 
0.5250 E 01* 
0.21*00 E 01* 

0.7800 E 03 
0.3500 E 03 
0.1900 E 03 
0.1108 E 03 
0.5383 E 02 

0.1*083 E 02 
0.2958 E 03 
0.2238 E 03 
0.7580 E 02 
0.3281* E 02 

0.101*1* E 02 
0.1*1*22 E 01 
0.2268 E 01 
0.1351* E 01 
0.5826 E 00 

0.3083 E 00 
0.1083 E 00 
0.6055 E-01 


TABLE 1. WEIGHTED ATTENUATION AND ABSORPTION COEFFICIENTS 
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